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1 Introduction 

The spontaneous breaking of chiral symmetry in Quantum Chromodynamics (QCD) is 
an important non-perturbative feature that determines the low-energy behaviour of the 
theory in a crucial way. In chiral perturbation theory, an effective description of QCD at 
low energies, the effect of spontaneous chiral symmetry breaking is encoded in the quark 
condensate £ and the Goldstone boson fields as the relevant degrees of freedom [1-3]. 
The effective theory then describes their dynamics while being strongly constrained by the 
requirements of chiral symmetry. In a specific regime of the theory, the so-called e-regime 
[4-7] where the mass of the Goldstone pion is small compared to the linear extent L 
of the space-time volume V = L 4 , i.e. m^L <C 1, the finite volume partition function is 
dominated by the zero-momentum mode of the pion field. In this situation, the topology of 
the gauge fields becomes important and as a consequence the partition function naturally 
separates into sectors of fixed topology [8] . 

Intriguingly, in the e-regime there exists another effective description of QCD which 
is based on random matrices [9-11]. In Random Matrix Theory (RMT) one models the 
fluctuations and correlations of the energy levels of the Hamilton operator of a physical 
system by considering the eigenvalues of matrices with random entries, but obeying the 
same global symmetries. For describing the Euclidean QCD partition function, for example, 
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one may construct a chiral random matrix model for the massless QCD Dirac operator. It 
turns out that in the microscopic regime this leads to the proper accumulation of eigenvalues 
at the origin, a characteristic feature of spontaneous chiral symmetry breaking. By now, 
one has a complete understanding of how and why chiral RMT represents an exact limit of 
the light quark partition function of QCD in the e-regime and a precise mapping between 
QCD, standard chiral perturbation theory and chiral RMT is established (see e.g. [12, 13] 
for a review and further references). As a consequence, RMT provides an interesting 
alternative approach for studying the universal properties of low-energy QCD. 

What is common to both effective theories is the fact that their intrinsic parameters, the 
effective low-energy constants, need to be determined externally, either from phenomenol- 
ogy, or directly from QCD using calculations from first principles, e.g. lattice QCD. In 
the past, there have been several efforts to extract physical parameters from numerical 
simulations in the e-regime, but all of them have been restricted to the use of lattice Dirac 
operators which obey an exact chiral symmetry [14-16]. Recently it has become evident, 
however, that simulations in the e-regime need not exclusively be reserved to these actions, 
but can also be done using standard Wilson fermions for which the chiral symmetry is 
broken by the discretisation. It has been understood in the context of Wilson chiral per- 
turbation theory (WChiPT) how to include the effects of the Wilson lattice discretisation 
in the standard regime of the theory [17-19] and this has also been achieved in the e-regime 
[20, 21]. Remarkably, also chiral RMT can be extended to include the effects of the Wilson 
lattice discretisation, leading to Wilson chiral Random Matrix Theory (WRMT) [22, 23]. 

Subsequently, a series of interesting statements have been made about the spectral 
density of the Wilson Dirac operator using WChiPT and/or WRMT for both Nf = 
(quenched) and Nf = 1 QCD [22-28]. One motivation is that a good description of the 
Wilson Dirac operator spectrum is most important for lattice simulations towards the 
chiral limit at fixed lattice spacing. In order to perform such simulations, a thorough 
understanding of the quark mass, volume and lattice spacing dependence of the spectrum 
is crucial for avoiding potential numerical instabilities. In that respect, the quenched 
approximation which we are studying in this paper can still be useful for Nf > simulations 
of QCD where m'EV 3> 1, because in that case the spectrum of the massless Dirac operator 
behaves essentially as in quenched QCD. 

Since topology plays a crucial role in the e-regime it is important to have good control 
over the topological properties of the system. A natural definition for the topological charge 
of the gauge fields is provided by the Atiyah-Singer index theorem which can be carried 
over to the lattice for chirally symmetric lattice Dirac operators [29]. One problem with 
simulating Wilson fermions in this extreme regime, though, is related to the fact that the 
chiral index is not well defined for the Wilson Dirac operator. 

However, there exists an exact relation between the number of real modes of the 
Wilson Dirac operator and the index as defined from the overlap operator [30, 31] using 
the Wilson Dirac operator as its kernel. One way to define a topological index for the 
Wilson Dirac operator is then provided by the eigenvalue flow of the hermitian Wilson Dirac 
operator, the so-called Wilson flow method [31-35]. The definition becomes unambiguous 
in the continuum limit and may hence serve as a natural definition of the index. It is an 
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astonishing result of the WRMT descriptions [22-24] that they can actually make definite 
statements about the hermitian eigenvalue flow and therefore make a natural connection 
to the various topological charge sectors of QCD. 

In this paper we investigate the spectral properties of the Wilson Dirac operator in 
quenched QCD in the e-regime corresponding to the microscopic regime of WRMT by 
means of numerical simulations. We calculate the low-lying spectrum of the hermitian 
Wilson Dirac operator and compare the distributions of the single eigenvalues to the ones 
from WRMT. Once the matching between the QCD data and WRMT is accomplished, we 
can check the predictive power of WRMT by scaling and comparing the results at different 
masses and volumes. For distinguishing the topological charge sectors we determine the 
Wilson flow on each gauge field configuration. As a by-product, the Wilson flow also yields 
the distribution of the real modes of the Wilson Dirac operator and this can again be 
compared to the chirality distributions as described by WRMT. 

2 Wilson random matrix theory 

We start by briefly recalling the WRMT that is expected to describe the spectrum of 
the Wilson Dirac operator. The Symanzik expansion has three operators - Wq, W7 and 
Ws - that provide the leading order corrections for WChiPT. Arguments from a large- iV c 
expansion indicate that the terms proportional to Wq and W7 are suppressed [36] and these 
terms are often set to zero in a first approximation. 
We consider matrices of the form 

W yiW^ a 8 B J ' V ; 

where A = and B = B* are (n + v) x (re + v) and re x n square matrices, respectively, 
while W is an arbitrary complex (n-\-v) x re rectangular matrix [23, 37]. We use tildes to 
distinguish the quantities from the analogous ones in the quantum field theory setup. The 
block structure of can be interpreted as a chiral decomposition with additional terms 
oc ds mixing the chiral sectors and thereby incorporating the effects of the operator Ws- 
Furthermore, the structure in eq.(2.1) guarantees that there are at least \u\ real eigenvalues. 
We are interested in the eigenvalues of the hermitian matrix 

DUrh) = %(D^ + rh) (2.2) 

where 75 = diag(l, . . . , 1, —1, . . . , —1) with (re + u) entries of 1 and n entries of —1 along 
the diagonal. 

The matrix elements of are randomly distributed according to the Gaussian weight 

P(A, B, W) = e -f Tr[^ + ^]-f Tr[Wt] ^ 

with N = 2n + v. The partition function Z^ for the WRMT is then given by integrating 
the matrices A, B, W over the complex Haar measure with the weight P(A, B, W), i.e., 

Z V N f (m,z;a 8 ) = I VAVBVWdet(D^ + m + z%) N f P(A,B,W) . (2.4) 
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The additional operators introduced by the Wilson lattice discretisation can be incor- 
porated by extending the partition function using the definition 

/oo y?i yj 

dy 6 dy 7 e ^ Z u Nf (m - y 6 , z - y 7 ; d 8 ) . (2.5) 
-oo 

Finally, the microscopic limit of WRMT is reached in the limit N — > oo while 

rhN, zN, diN 1/2 (2.6) 

are kept fixed. It is therefore natural to think of ./V as corresponding to the space-time 
volume in the quantum field theory setup. One can show that the WRMT partition function 
reproduces the leading terms in the e-expansion of WChiPT [23] and the corresponding 
low-energy constants can be identified as 

mN otmZV, 2N oc zZV, a 2 N oc a 2 WiV. (2.7) 

Note that the WRMT setup described here can only reproduce WChiPT with a fixed 
choice for the signs of the low-energy constants. Using the conventions of [23], ds will map 
into a positive W%, while non-zero values a§ and d 7 imply negative values for W% and Wj. 
From the corresponding partition functions analytic expressions for the distributions of the 
eigenvalues of the QCD Dirac operator in the e-regime can be obtained both in WChiPT 
[23] and in WRMT [24]. 



3 Lattice setup 

Before discussing our results we also need to specify the setup for our lattice calculation. 
The massive Wilson-Dirac operator can be written as 

D w {m)= l - lv {V v + Vl)- l -VlV u + m, (3.1) 

where V„, V* denote the covariant forward and backward lattice derivatives, respectively, 
7jy are the Euclidean Dirac matrices and m is the bare quark mass. The Wilson-Dirac 
operator is 75-hermitian, 

D\, = 75J Dw75 (3.2) 

and hence 

Ds = 75£>w("i) ( 3 - 3 ) 

is hermitian, = D§. The dependence of the eigenvalues A|(m) of D§(m) on m defines 
the eigenvalue flow and from first order perturbation theory one can infer 

The eigenvalues of Z?\y are denoted by A^. Since the Wilson Dirac operator is non- normal, 
i.e. [D^,!)^] 7^ 0, one needs to distinguish the left and right eigenvectors which define a 
bi-orthogonal system. As a consequence the eigensystems of Z?w and -D5 are related in a 
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non-trivial way, except for the subset of Ds-eigenvalues which are exactly zero for a given 
value of m. In that case one has 



75(Z?w + m )ip = <^=^ Dy/ip = —mip , (3.5) 

i.e., a real mode £ K of -D-yv corresponds to a zero mode of D^{m = — A^ ) and the 
chirality of the mode is given by the derivative of the Wilson flow, eq.(3.4), at m = — A^. 

Here we are specifically interested in the chirality distribution of the real modes A^ of 
Dyy defined as 

PxW = ( E 5 ^ + A)sign«% 5 |fc» ) (3.6) 



and the distribution of the k-th eigenvalue A|(m) of D§ given by 

pl(X,m) = (5(XUm)-X)) . (3.7) 

These are the distributions which we compare to the analytic expressions obtained from 
WRMT. Up to now only the expression for the full distribution p^ = J2k p\ 1S known, 
so the single eigenvalue distributions p\ are obtained from WRMT directly by numerical 
simulation. 

3.1 Definition of the topological charge 

A crucial ingredient for the proper treatment of the configurations in the analysis is their 
correct assignment to the various topological sectors. As discussed in [23] summing up the 
signs of the chiralities of the real modes A^ € M of Z?w f° r a given gauge field configuration 
yields an index of the Wilson Dirac operator, 

v= E sign«% 5 |fc». (3.8) 

Note that the sum over the real modes is actually only a sum over the real modes of 
Dw(m) in the neighbourhood of m, i.e. it excludes the doubler real modes. Since the real 
modes are in correspondence with the zero modes of D§, this index is in agreement with 
the topological charge as defined through the Wilson flow method [31-35]. In practise, the 
Wilson flow yields all the real modes of Z?w up to a cut-off m cut , so the practical realisation 
of eq.(3.8) is given by 

u= Yl sign«fe| 75 |fc», (3.9) 

Ajf <m cut 

where the restriction takes into account only the real modes up to a certain cut-off m cu t- 
Consequently all real modes with A]^ > m cut are considered to correspond to doubler 
modes. In that sense the topological charge is not well defined at finite lattice spacing - as 
it should be - but becomes so in the continuum limit. This ambiguity is not a surprise and 
is reflected in any definition of the topological charge at finite lattice spacing. In particular, 
also the definition of the topological charge via the index of the overlap Dirac operator [30], 

D = 1 + 75 sign[7 5J D K (-m)] , (3.10) 
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is ambiguous, since it depends on the mass shift — m entering the definition of the hermitian 
kernel Dirac operator 75-Dk- The correspondence between the index of the overlap operator 
and the charge from the eigenvalue flow is exact, so m in the overlap operator corresponds 
to m cut in the definition eq.(3.9). 

While the chirality of the real modes is ±1 in the continuum, for the non-normal Wil- 
son Dirac operator this holds only approximately (cf. also [38, 39]). The non-normality of 
Z?w is a consequence of its specific symmetry structure and can not be avoided, however, 
one can reduce it, e.g., by modifying the covariant derivatives. In our study we use co- 
variant derivatives with various levels of HYP-smeared gauge fields, since this is a simple 
but efficient mean of suppressing the non- normality [40]. The effect of the smearing is 
manifold. Most importantly, it improves the separation of the physical zero modes from 
the doubler zero modes and in this way reduces the ambiguity of the charge definition, 
i.e. its dependence on m cut . As a side effect, the chirality of the zero modes is improved 
significantly with increasing smearing level. 

Obviously, any remaining ambiguity in the topological charge assignment is due to ad- 
ditional crossings beyond the value of m cu t employed in the charge definition. Interestingly, 
WRMT can make predictions about the number of these additional crossings [24] and we 
will investigate this in section 4.3.1. 

Based on these considerations, we present in the following the results for a Wilson Dirac 
operator with covariant derivatives based on gauge fields with 5 HYP-smearing levels, and 
use eq.(3.9) as our definition of the topological charge with m cu t as given in table 1. 

3.2 Simulation parameters 

We have simulated several lattices with various sizes and lattice spacings in Nf = QCD 
using the Wilson gauge action. The simulation parameters for these quenched ensembles 
are described in table 1. Within the group A we have kept the lattice spacing constant in 
order to study potential finite volume effects. In contrast, for ensembles A\ and B\ the 
volume V = L 3 x Lt is roughly kept constant in physical units. Table 1 also provides a 
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0.50 


0.12(13) 



Table 1: Simulation parameters for the Nf = ensembles together with the value of m cut used in the 
definition of the topological charge. 

numerical estimate of the derivative of the charge with respect to the cut employed in the 
Wilson flow. A graphical presentation of this quantity as a function of the mass is shown 
for each of the ensembles in figure 1 . The derivative in each of the ensembles shows a sharp 
peak in the neighbourhood of the critical mass, after which a more or less rapid decay sets 
in depending on the volume and the lattice spacing. In absolute terms, the derivative of 
the charge is practically negligible at the point where the mass cut is introduced. Using 
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Figure 1: Derivative of the absolute charge with respect to m cut for the ensembles described in table 1. 

a linear extrapolation through the tail of the distributions as a rough estimate, we should 
expect to be misassigning up to half a percent of the configurations in ensemble B\ and at 
most about one percent for the ensembles A\ and A^- 

4 Results 

4.1 Cumulative hermitian eigenvalue distributions 

We start by considering the distributions /05(A) of the eigenvalues of the hermitian Wilson 
Dirac operator D§ and compare them with the analytic expressions. So far the expres- 
sions from WRMT and WChiPT comprise only the full distributions for the lowest lying 
eigenvalues. In figure 2 we show the full distributions for our ensemble A\ in the charge 
sectors v = and v = 1, together with fits to both that will be discussed in more detail 
below. Some of the expected features of the spectrum are apparent, the most prominent 
being the enhanced peak on the positive side for v = 1. However, it is clear - also from 
a comparison of the data to the relatively delicate features of the fitted lines - that it 
would be rather challenging to push the analysis beyond a qualitative comparison on the 
basis of these distributions alone. Luckily, our data obtained from the simulations contain 
much more information than just the full distribution, since we have access to each single 
low- lying eigenvalue. Indeed it turns out that the shapes of the single eigenvalue distri- 
butions are rather sensitive to the low-energy constants Wq, Wj and Wg. For this reason, 
our fits concern single eigenvalue distributions obtained from explicit WRMT simulations 
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Figure 2: The density of the lowest lying eigenvalues in the topological charge sectors v = (top) and 
v = 1 (bottom), for the ensembles Ai. Blue lines are fits accounting just for Wa, red lines are fits including 
the effects of W% and Wr. These lines were calculated by means of WRMT, using 250k samples. 

including the effects of W% and possibly Wq,Wi. The fit parameters in WRMT are then 
rh, og, 07, as and YiV which sets the scale and connects the WRMT to the QCD data. 

Binning the data into histograms introduces an unwelcome dependence on the bin size 
and hence we choose to compare cumulative distributions. In order to assess the quality of 
the fits we employ the Kolmogorov-Smirnov (KS) test statistics. If Sqcd{x) and /Swrmt^) 
are the estimates of the cumulative distribution functions obtained from the QCD and the 
WRMT simulation, respectively, the KS statistics is given by 

d max = max |Sqcd(z) - Swrmt(z)| • (4.1) 

— oo<a;<oo 
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Since we fit several distributions at once, it is not clear how to translate <i max into a quality 
of fit measure and we will restrict ourselves to reporting only the value of d max - 

Finally, the errors on the distributions we present in the following are obtained by 
a bootstrap procedure and for better visibility represent a 95% confidence bound, rather 
than a regular standard deviation. 

A summary of all the fit results discussed in this paper is given in table 2 for convenient 
reference. 
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Table 2: Fit results for the fits of WRMT to the QCD ensembles Ai, A2 and B\ in the various topological 
sectors \v\. N con { denotes the number of configurations used in the fits. The bare mass m of the Dirac 
operator and T,V are given in lattice units, while the remaining columns are the fitted dimensionless WRMT 
parameters. The total number of the lowest-lying eigenvalues fitted on both the positive and negative side 
of the spectrum is iV cv = 4 for fit F 1 and iV cv = 8 for all other fits. 



4.1.1 Influence of We and W7 

As mentioned above, large- N c arguments indicate that the terms proportional to Wq and 
W7 are suppressed [36]. In order to assess the importance of these low-energy constants 
we compare fits including only W$ (fits F5 and F§ in table 2) to ones which include all 
three low-energy parameters Wq,Wj and Wg (fits F 2 and F%). The resulting distributions 
obtained in ensemble Ax in charge sector v = at m = —0.07 are displayed in figure 3 and 
a visual examination makes it clear that the fits including Wq , W7 are superior in describing 
the shapes of the single eigenvalue distributions. This is corroborated by the values of d max 
in table 2 which grow by about 30-50% when Wq,Wj are excluded. So while W% appears 
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to account for large parts of the deviation from the chiral setup, we do find Wq and W7 to 
be important as well and henceforth include them in all further fits. 
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Figure 3: Cumulative eigenvalue distributions for ensemble Ai in charge sector v = compared to WRMT 
fits including contributions from low-energy constants We, W7, Wg (fit F 2 , upper plot) and only Wg (fit F5, 
lower plot). Blue lines denote fitted distributions from WRMT, red lines are the predictions following from 
the fit. 

Fits including and excluding the effects of Wq and W7 have also been added to figure 2. 
Here one can clearly observe the influence of the additional operators, that work to largely 
wash out the characteristic but subtle undulating pattern observed on top of the spectrum 
as predicted by the WChiPT expressions [23]. Our observation that fitting these added 
parameters tends to produce non-vanishing values implies that one should not necessarily 
expect to observe such a pattern in practise. This once again underlines the importance of 
examining the distributions of the separate eigenvalues rather than the overall distribution, 
as the apparent impact of these operators would imply few distinct features of the spectrum 
remain to constrain fits by. 

4.1.2 Sensitivity to topology and predictive power 

Next we investigate the sensitivity of the WRMT fits to the various topological sectors 
and the predictive power of the fits as a first non-trivial test. A possible way to proceed is 
to fit the distributions of a selected set of eigenvalues in a given topological sector to the 
distributions obtained from WRMT and check the consistency for the other eigenvalues 
and in the other topological sectors. 

In figure 4 we show the results of such an exercise for the spectrum of -D5 at a bare 
mass m = —0.07. We choose to fit the distributions of the lowest two positive and negative 
eigenvalues in the charge sector v = 1, since this is where we have the largest statistics. 
As before, the fitted cumulative WRMT distributions are displayed as the blue lines, while 
the red lines are the predictions for the remaining eigenvalue distributions following from 
the fit. The results for the fit parameters and the value of d max can be found in table 2 in 
the row labelled by F\. 
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Figure 4: Cumulative eigenvalue distributions for ensemble A\ in charge sectors v = 0, 1 and 2 from top 
to bottom. The QCD data is displayed in black. Blue lines denote the fitted distributions from WRMT 
(fit F\ in table 2), while red lines are the predictions following from the fit. 



A number of observations can be made at this point: 

(1) The four eigenvalue distributions in the charge sector v = 1 can be described very well 
by WRMT, both in terms of the position and the shape, as reflected by the low value of 
dmax reported in table 2. 

(2) The distributions in the charge v = sector are all predictions. Again, the agreement 
is very good up to about |A|SV ~ 12 where deviations in the tail of the distributions away 
from start to show up. One may assign this to the onset of the bulk which could be 
related to the Thouless energy where the RMT description is expected to break down. 

(3) For the predicted distributions in the charge v = 1 sector we find that WRMT tends 
to underestimate the tails of the eigenvalue distributions away from 0. 

(4) The same effect is visible in the other predicted eigenvalue distributions in charge sector 
v = 2, but in addition we now also find a mismatch in the position of the distributions. 

(5) It seems that both these mismatches are more severe on that side of the spectrum 
where the would-be real modes are accumulating. One may conclude from this that lattice 
artefacts seem to be pushed from the would-be real modes to the nearby modes in the 
spectrum, and that this effect is enhanced for increasing topology. 
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Figure 5: Cumulative eigenvalue distributions for ensemble A\ with volume Vi = (1.52fm) 4 together with 
the best fits to WRMT. Charge sector v = (top, fit F 2 ) and v = f (bottom, fit F :i ). 

4.1.3 Scaling of the distributions 

An important prediction of the WRMT is the universal scaling of the distributions with 
the volume, as implied by equation 2.7, and with the mass. 

We first investigate the scaling of the distributions with the volume. For this purpose 
we compare our best fits -F2-3 on ensemble A\ with a volume V\ = (1.52fm) in the charge 
sectors v = 0, 1 to the corresponding fits -F9-10 on ensemble A2 with volume V2 = (1.28fm) 4 . 
Both data sets are obtained at m = —0.07. The fits are displayed for visual inspection in 
figure 5 for A\ and figure 6 for A2. We note that the quality of the fits on the smaller 
volume is clearly deteriorating, with d max -values increasing by a factor 3 in charge sector 
v = and 4 in sector v = 1. We assign this to the fact that the volume V2 = (1.28fm) 4 
is clearly too small for the system to be in the e-regime. Nevertheless, with the results in 
table 2 we may quantitatively check the scaling of the parameters. We expect 0,5,0,7 and 
d 8 to scale like (V1/V2) 1 / 2 = 1.41 and T.V like V1/V2 = 2.0. Except for a 8 this is what we 
roughly observe at least in the charge sector v = where the small volume is indeed better 
described by WRMT. 

Next we investigate the scaling of the distributions with the mass. For this purpose we 
plot in figure 7 the distributions obtained on A\ in the charge sectors v = 0, 1 at a different 
value of the bare mass of m = —0.055. These fits should again be compared with those 
in figure 5. From a comparison of the fit parameter rh, we can infer a value of the critical 
mass m cr it ~ —0.09. As we will see in section 4.3 this matches rather well with what we 
find from the distributions of the real modes. On the other hand, we observe a rather poor 
quality of the fits Fj and F% as compared to F2 and F3. Of course, the increase in mass 
will push all the eigenvalues in the spectrum towards the bulk making it more difficult for 
the WRMT to describe them properly. Indeed, this is also reflected in the enhanced errors 
on all the fit parameters impeding a more quantitative comparison of the two spectra at 
different masses. 
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4.2 Eigenvalue averages 

With respect to the full eigenvalue distributions discussed in the previous section, averages 
of single eigenvalues represent a strong reduction of the amount of information available 
for a comparison. Nevertheless, since the averages can be determined to high statistical 
precision, a comparison between the QCD data and the WRMT description at this level 
may indicate systematic deviations more clearly. Figure 8 summarises the results of our 
fits to the two main ensembles A± and B\ in terms of these averages. The plots show the 
averages for each separate eigenvalue, both as measured on the lattice and as produced by 
the fitted WRMT. In each topological charge sector, values on the left denote the negative 
eigenvalues, while the values on the right give the positive ones. For the charge sector 
v = the two sides are symmetric and have been merged. 

We find that the WRMT appears to work very well for trivial topology, both for the 
coarse ensemble B\ and the smooth one A%, as can also be inferred from table 2. For higher 
charges, however, the eigenvalues with chiralities opposite to the sign of the topological 
charge are systematically pulled towards zero, while the ones with the same chirality sign 
come out higher than what WRMT tends to predict. This points to the mixing between 
eigenvalue sectors of different chirality being stronger, and hence lattice artefacts being 
larger than what WRMT is able to describe. The effect seems to be more pronounced in 
the higher charge sectors and at coarser lattice spacing. This is corroborated by the fact 
that for the ensemble Ai the charge sector v = 1 is still well described, whereas for B\ 
large deviations start to appear already there. 

4.3 Chirality distributions 

The determination of the topological charge through the Wilson flow provides us with all the 
zeros of D^(—m) up to the cut-off m cut together with their associated chirality. According 
to the relationship in eq.(3.5), the locations of these zeros can be directly related to the 
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Figure 6: Cumulative eigenvalue distributions for ensemble A2 with V2 = (1.28fm) 4 together with the 
best fits to WRMT. Charge sector v = (top, fit Fg) and v = f (bottom, fit F 10 ). 
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Figure 7: Cumulative eigenvalue distributions for ensemble Ai at bare mass m 
v = (top, fit Ft) and v = 1 (bottom, fit Fs). 
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Figure 8: Comparison of the measured and fitted eigenvalue averages, for the ensembles A\ and B\ in the 
different topological charge sectors. Negative eigenvalues are plotted to the left, positive eigenvalues to the 
right. Black points indicate the measured data, while red lines indicate the values from WRMT based on 
the fits F 2 , F 3 , F 4 , Fu and F 12 from table 2. 



real eigenvalues of the operator Z?w- The latter can also be directly obtained from 
WRMT, once the parameters have been determined, allowing for the comparison of the 
chirality distributions. 

Figure 9 shows a comparison of the chirality distributions that were measured for our 
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ensembles A\ (in red) and B\ (in blue). Empty bars show the number of eigenvalues 
with positive and negative chirality, while the filled bars give the net chirality distribution. 
To simplify the comparison, the measured values of A = — A^ were multiplied by EV as 
determined from our fits and shifted such that the peak of the distribution is located at 
in a rough approximation of the critical mass. The distributions for B\ were rescaled to 
compensate for differences in the statistics per charge sector. 

The distributions for both ensembles A\ and B\ are sharply peaked around m c , but 
long tails remain towards more negative values of A. It is clear that the distribution is more 
sharply peaked, and the tail shorter, for the ensemble A± with its significantly smaller lattice 
spacing. In addition, the number of cancellations is significantly higher for the ensemble 
B\. The occurrence of conjugate pairs of real modes with opposite chirality, responsible 
for these cancellations, will be suppressed towards the continuum limit. This comparison 
therefore shows a clear decrease in discretisation artifacts. 

Analytical results for the chiral density [23] based on WChiPT tend to show a distinct 
structure of merged peaks, the number of which is given by the topological charge. This is 
in contrast to what we observe in any of our ensembles. One possible cause of this could be 
the influence of the Wq and W7 operators, as already indicated in figure 2. To investigate 
this further, we obtained the real eigenvalues from 100k samples of the WRMT with v = 2. 
We used sets of parameters matching our A\ ensemble, both including clq and ay. The 
chirality distributions obtained from this calculation are provided in figure 10. The effect 
of the additional operators is clearly observable, as the peak is broadened and all internal 
structure wiped out. Given the non- vanishing values we find for the parameters a,Q and 07 
in our fits, we should not expect to see any structure in our chirality distributions. 

While real eigenmodes can and do appear for the topologically trivial sector, a chirality 
distribution equal to zero is guaranteed because of the symmetry between positive and 
negative values of the topological charge. In order to facilitate a quantitative comparison 
between the chirality distributions of QCD and the ones from WRMT in the sectors with 
non-trivial topology we again consider the cumulative distributions. Figure 11 shows those 
found in ensemble A± from configurations with v = 1 (red) and 2 (blue) in solid lines, with 
the WRMT predictions in dashed lines. Note that in this plot, m cut would correspond to 
a value of A = —0.3. 

We find a clear distinction between the QCD and the WRMT distributions, as the 
QCD results are asymmetric around 0. It is interesting, however, to observe that the 
distributions do in fact match quite well up to the determined value of the critical mass. 
Beyond that point, the chirality distribution as described by WRMT drops to zero quite 
rapidly, while crossings keep occurring for the QCD theory. The fact that the derivative of 
the absolute charge with respect to the mass cut-off has become very small for larger values 
of the cut-off (cf. table 1) shows that this deviation cannot be attributed straightforwardly 
to a lack of statistics or misassigned configurations. On the other hand, the structure of 
the WRMT prevents any of the operators included from breaking the symmetry of the 
distribution around 0, so the mismatch must be due to lattice artefacts in the chirality 
distributions which appear not to be described by WRMT. 
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Figure 9: Chirality distributions for the three lowest charge sectors of ensembles Ai (red) and Bi (blue). 
The distributions have been shifted for convenient comparison and have been normalised to contain the same 
statistics per charge sector. The dotted line corresponds to the value of m cu t employed for the calculation 
of the charge in ensemble Ai. 



4.3.1 Additional modes 

Apart from the chirality distribution, WRMT and WChiPT allow for inferring a relation- 
ship between the charge sector and the average number of additional crossings N^dd, i- e -> 
the difference between the total number of crossings found in the Wilson flow and the 
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Figure 10: Chirality density as calculated by WRMT for v = 2, excluding (red) and including (blue) the 
effect of as and 117. 
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Figure 11: Cumulative distribution function for the chirality distribution of the A\ ensemble. Red lines 
refer to v = 1, blue lines to v = 2. Solid lines are obtained from the QCD simulation, while dashed lines 
indicate the expectation from WRMT. 



absolute charge of the configuration. For small lattice spacing a, this relationship is given 
in [24] as 

iV add oca 2 ^ +1 ), (4.2) 

a pure power law. The determination of iV a dd becomes progressively harder for larger val- 
ues of the charge, as both the number of such configurations and the expectation value 
of -/V a( id decrease rapidly. Nevertheless, the data for ensembles Al and Bl displayed in 
figure 12 appears consistent with the predicted power-law behaviour. Moreover, the ratio 
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of the coefficients of the linear terms should be given by the ratio of the lattice spacings. 
Performing linear fits to the initial three (for A\) and four (for B{) points, we obtain the 
coefficients -0.785 and -0.471, respectively. The ratio then comes out as 1.67, to be com- 
pared to 1.64 from the ratio of lattice spacings as given by ro/a. This excellent agreement 
may be somewhat accidental here, given the limited number of reliable points and the clear 
sensitivity to the fit range. Nevertheless, we conclude that there are strong indications that 
the power law for N^dd derived in [24] holds for our data, in spite of the limited agreement 
with the prediction for the chirality distribution and the somewhat coarse lattice spacing 
in B x . 



T3 1 



bJO 
O 




Figure 12: The logarithm of the average number of additional crossings iV a dd versus the absolute charge 
\v\. Shown are values for ensemble A\ (red) and B\ (blue). 



5 Conclusions 

We have investigated the spectral properties of the Wilson Dirac operator on the lattice for 
Nf = quenched QCD in the e-regime by numerical simulations. The resulting eigenvalue 
distributions /05(A) of the hermitian Wilson Dirac operator D$ and the chirality distribu- 
tions p x (X) of the real modes of the Dirac operator L>w can be compared to expressions 
resulting from Wilson chiral perturbation theory (WChiPT) in the e-regime, or equiva- 
lently by Wilson random matrix theory (WRMT) in the microscopic limit. In order to 
do so, it is important to have good control over the different topological sectors of QCD 
and we have taken care of this by employing the topological index based on the eigenvalue 
flow of -D5 as the definition of the topological charge. Towards the continuum limit this 
definition becomes equivalent to the index of the overlap operator. 

Concerning the eigenvalue distributions /55(A) of the fc-th eigenvalue, we find that they 
can be described rather well by the corresponding distributions as given by the WRMT. It 
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turns out that by fitting these distributions, it is indeed possible to extract the leading low 
energy constants Wq , Wj and Ws for Wilson fermions from the eigenvalue spectra of the 
Wilson Dirac operator in the e-regime, on volumes as small as V ~ (1.5fm) 4 and lattice 
spacings as coarse as a ~ O.lfm. However, in order to describe the shape of the distributions 
properly, it is in fact necessary to include the effects proportional to Wq and W7 in addition 
to the ones proportional to Ws- From the resulting fits, one can extract values for Y>V at 
the one percent level. However, since these values should rather be regarded as effective 
parameters to match the spectral distributions of D§ with those of -D5, we refrain from 
quoting a physical value for the corresponding quark condensate. 

Our results indicate that the description of the microscopic Wilson Dirac operator 
spectrum by WRMT works best in the topological charge sectors with low values of the 
charge and starts to break down as the topological charge increases, the more so the smaller 
the physical volume. In addition, the chirality distributions obtained from our simulations 
show an asymmetry around with an enhanced tail for > which can not be described 
by the operators included in WRMT and WChiPT parametrising the leading order lattice 
discretisation effects. It appears that the mixing of eigenmodes of different chiralities due 
to the breaking of chiral symmetry by the Wilson term is stronger than anticipated by the 
description of WRMT and WChiPT, at least at the lattice spacings investigated in this 
work. Nevertheless, an important result of WRMT is that the number of additional real 
modes in excess of the topological charge v is suppressed for large values of v [24] and we 
do find that this is indeed the case. 

It would be particularly interesting to see how WRMT works for the special case of 
Nf = 1 QCD. In that situation there is no spontaneous breaking of chiral symmetry, 
but a WRMT description is nevertheless available. Since the theory defined with Wilson 
fermions suffers from a sign problem towards the e-regime, any information on the spectral 
properties would be useful in order to avoid numerical instabilities of such simulations. 
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